Local and nonlocal pressure Hessian effects in real and synthetic fluid turbulence 
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The Lagrangian dynamics of the velocity gradient tensor A in isotropic and homogeneous tur- 
bulence depend on the joint action of the self-streching term and the pressure Hessian. Existing 
closures for pressure effects in terms of A are unable to reproduce one important statistical role 
played by the anisotropic part of the pressure Hessian, namely the redistribution of the probabilities 
towards enstrophy production dominated regions. As a step towards elucidating the required prop- 
erties of closures, we study several synthetic velocity fields and how well they reproduce anisotropic 
pressure effects. It is found that synthetic (i) Gaussian, (ii) Multifractal and (iii) Minimal Turnover 
Lagrangian Map (MTLM) incompressible velocity fields reproduce many features of real pressure 
fields that are obtained from numerical simulations of the Navier Stokes equations, including the 
redistribution towards enstrophy-production regions. The synthetic fields include both spatially 
local, and nonlocal, anisotropic pressure effects. However, we show that the local effects appear to 
be the most important ones: by assuming that the pressure Hessian is local in space, an expression 
in terms of the Hessian of the second invariant Q of the velocity gradient tensor can be obtained. 
This term is found to be well correlated with the true pressure Hessian both in terms of eigenvalue 
magnitudes and eigenvector alignments. 

PACS numbers: 02.50. Fz, 47.53.+n, 47.27.Gs 



I. INTRODUCTION 

The study of the velocity gradient tensor in fully de- 
veloped turbulence has lead to interesting findings and 
has contributed to improved understanding of many sta- 
tistical and geometrical properties of turbulent flows. In 
particular, recent progress has been made in the study 
of the Lagrangian dynamics and modeling of the veloc- 
ity gradient tensor (see Ref. [l[ for an overview of the 
subject). This tensor is given by Ay = dui/dxj, where 
u is the velocity vector. Taking a spatial gradient of the 
Navicr-Stokes equations, the following transport equa- 
tion for A is obtained: 



— = A 2 P 

dt 



/AA 



(1) 



where d/dt stands for the Lagrangian time derivative, 
v is the kinematic viscosity and = d 2 p/dxidxj is 
the pressure Hessian. The first term —A is the self- 
stretching term. The restricted Euler (RE) approxi- 
mation, which assumes an isotropic pressure Hessian 
Pij = — tr(A 2 )5ij /3 and neglects viscous effects, leads to 
an autonomous set of coupled ordinary differential equa- 
tions Q . The intrinsic dynamics of the RE system leads 
to a finite time divergence of the components of A dur- 
ing which the vorticity u> = V A u gets aligned with 
the eigenvector of the rate of strain, S = (A + A T )/2, 
associated with the intermediate eigenvalue, as often ob- 



served in real turbulence 

0, S H- To prevent the devel- 
opment of unphysical finite time singularities, both the 
anisotropic part of the pressure Hessian and the viscous 
diffusion term have to be modeled. This was the subject 
of former works @-Q|- In particular, closures were pro- 
posed in Ref. Q for P and MA in terms of the local 
value of A. The local closures of Ref. 0, when inserted 
into the dynamics generated by Eq. (JTJ under the action 
of a stochastic forcing term, lead to stationary statis- 
tics of A along Lagrangian trajectories which compare 
well with those obtained from direct numerical simula- 
tions (DNS) of the Navier-Stokes equations at moderate 
Reynolds numbers ||. At higher Reynolds numbers, pre- 
dictions of the stochastic model proposed in Refs. @, Q 
turn out to become unrealistic, mainly because of the 
weakness of the closure for the anisotropic part of the 
pressure Hessian. 

Indeed, the pressure Hessian is related to the spatial 
distribution of the velocity gradient using singular inte- 
gral operators [§-[l3[ according to: 



d 2 p 



dxidx 



_ = -tr(A 2 )^ - P.V. J fcy(x - y)tr(A 2 )(y)dy. 



(2) 

Above, the integral is understood as a Cauchy principal 
value (P.V.) and kij is the Hessian of the Green's function 
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for the Laplacian operator, namely 



fcij(x) = 



1 



|x| 2 <5, 



dxidxj 47r|x| 



47r|x| 5 



(3) 



One can see from Eq. ([2]) that only the isotropic part 
of the pressure Hessian is purely local (the first term on 
the right-hand side (RHS) of Eq. ©). All the nonlocal 
effects of pressure Hessian enter through the anisotropic 
part (or deviatoric part corresponding to the second term 
in the RHS of Eq. ([2])). Hence, in this view, the RE ap- 
proximation can be understood as the neglect of all the 
nonlocal effects implied by the incompressibility condi- 
tion (or pressure field). 

In order to quantify the precise action of pressure in 
numerical turbulent flows, it was proposed, in Ref. @, to 
study the probability current associated with pressure in 
the plane spanned by the two highly relevant invariants 
of A, R and Q (the so-called i?Q-plane). One of these 
invariants, defined as, 



Q = -^tr(A 2 



-\u>\ 2 - -tr(S 2 ) 
4 2 v ' 



(4) 



quantifies the net balance, or competition, between en- 
strophy and dissipation. The other important invariant, 
defined as 



J?=-^tr(A 3 ) = 



1 1 3 

--UiSijUJj - -tr(S ), (5) 



quantifies the competition between enstrophy production 
and strain skewness (i.e. dissipation production). As it 
will be recalled in the following, in terms of the velocity 
gradient evolution in statistically stationary turbulence, 
pressure has two important roles. First of all, pressure 
counteracts the development of the singularity implied by 
the self-stretching term. This feature is found to be well 
reproduced by existing closures 0. The other impor- 
tant pressure action is the redistribution of probabilities 
towards enstrophy production dominated regions (i.e. to- 
wards R < 0). This is not reproduced well by existing 
closures. As discussed in Refs. [HHH], a related deficiency 
of the closures is that they all predict that the pressure 
Hessian is proportional to Q. In the R — Q plane dynam- 
ics, this implies that when Q = 0, the effect of pressure 
Hessian also vanishes. For real turbulence, there is no 
such vanishing of pressure Hessian effects when Q = 

a. 

In this article, we investigate whether these particu- 
lar features of pressure (redistribution of probabilities 
towards enstrophy production, and non- vanishing action 
even when Q = 0) is inherent to true Navier- Stokes tur- 
bulence, or can also be observed in various approxima- 
tions, namely synthetic turbulent velocity fields. Var- 
ious types of synthetic fields are considered. The first 
type of synthetic field considered is Gaussian fields ob- 
tained by superposing random-phase Fourier modes with 
prescribed spectra. The second type is called 'multifrac- 
tal' [H|, and consists of a Gaussian field whose vortic- 
ity field is amplified by means of the 'Fluid Deformation 



Closure' and made consistent with multifractality's long- 
range correlations in physical space [l6j]. The third type 
of synthetic field is gen erated using the Lagrangian map- 
ping technique (l7l Il8|. It also relies on random-phase 
Gaussian fields but then applies a multi-scale deforma- 
tion of fluid particles using a simple Lagrangian mapping. 
For each of these synthetic velocity fields, a pressure field 
is obtained numerically by means of the pressure Poisson 
equation. As will be seen, unlike the local closures dis- 
cussed above, these synthetic fields reproduce many cor- 
rect features of the pressure Hessian. In particular, they 
will be shown to reproduce the redistribution of probabil- 
ity towards enstrophy production, as well as displaying 
non- vanishing action, even when (5 = 0. 

The second part of the paper studies to what degree 
spatial locality is important in determining these prop- 
erties of the anisotropic part of the pressure Hessian. As 
can be seen from the expression for the pressure Hessian 
(Eq. fl2J)), the anisotropic part of the pressure Hessian 
is also the part that is spatially nonlocal, i.e. the part 
that requires knowledge of tr(A 2 ) at positions y / x. 
Arguably, the more non-local effects are important, the 
more challenging it is to formulate closures in terms of lo- 
cal quantities. In order to examine the degree of locality, 
in the second part of this paper we decompose the space 
integration in Eq. ^ into two parts, a local part given 
by the integration over a small ball of radius given by 
the Kolmogorov length scale tjk, and the remainder be- 
ing the 'nonlocal' portion. We will show that neglecting 
the second non-local contribution leads to an expression 
that models the anisotropic part of the pressure Hessian 
in terms of the Hessian of the invariant Q. Using DNS 
data, this expression is compared with the true pressure 
Hessian. 



II. PRESSURE HESSIAN FROM DNS AND 
SYNTHETIC VELOCITY FIELDS 

A. Probability current in the i?Q-plane 

1. Definition of the probability current 

We follow the approach used in Refs. [1, based 
on a Fokkcr-Planck equation for the dynamics of R and 
Q. To summarize the approach, we remark that it can 
be shown that the time evolution along a Lagrangian 
trajectory of the non-dimensional invariants R* = R/a 3 
and Q* = Q/cr 2 is given by 

^ = -3R*-±A lk Hl-±A tk H» kt and (6) 



if = I (Q * )2 ~ h AikAklH ^ ~ h AikAklH ^ ' (7) 

where a 2 = (tr(S 2 )) is the strain variance, and t* = at is 
the non-dimensional time. Also, H p stands for (minus) 
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the deviatoric part of the pressure Hessian, i.e. 



/ d 2 p 



d 2 p 



\dxidxj 3 dxkdxk 



(8) 



and = vV 2 A. is the viscous term (recall that in the 
RE approximation, LP = H" = 0). The Fokkcr-Planck 
equation describing the time evolution of the joint density 
V(Q*,R*) may be written as: 



dV 

dt* 



d 

ao* i yy = , 

dR* 



(9) 



where the divergence of the probability current W con- 
trols the time variations of the joint probability density 
V . The probability current can be written in terms of 
conditional averages as 



W 



dQ* 
ell • 

.//,' 

dt* 



Q*,R*)V(Q*,R*) 



(10) 



and can be decomposed into W = Wre + W p + W v , 
where the probability currents Wre, W p and W v are 
associated, respectively, to the effects on the Lagrangian 
evolution of the invariants Q and R (Eqs. © and ([7])) 
of the restricted Euler term —A 2 , of (minus) the pres- 
sure Hessian — P and of diffusivity ^V 2 A entering in Eq. 
(TT]). In this article, we will focus on the probability cur- 
rent associated with the pressure Hessian W p - It can be 
written as 



W v 



-A lk H p kl /a 3 
-A lk A kl HfJo A 



Q*,R*)V(Q*,R*) . (11) 



More details are provided in Ref. 



2. DNS velocity fields 

In the following, we will make extensive use of data 
from standard direct numerical simulation (DNS) of the 
Navier-Stokes equations, for a Taylor-based Reynolds 
number of order 1Z\ = 145. DNS is based on a 
pseudo-spectral method with 2nd-order accurate Adams- 
Bashforth time stepping; the computation box is cubic 
(size 2n) with periodic boundary conditions in the three 
directions and spatial resolution 512 3 . Statistical station- 
arity is maintained by an isotropic external force acting 
at low wavenumbers in order to ensure a constant energy- 
power supply. It provides, in the units of the simulation, 
a constant energy injection rate (e) = 0.001. The kine- 
matic viscosity of the fluid is v = 0.000285. The Kol- 
mogorov's scale is r/K = (zy 3 /(e)) 1 / 4 = 0.0123 so that 
dx/rjK ~ 1, since dx = 2n/512. 

We display in Fig. QTa) the vector plot and stream- 
lines of the probability current W p associated with the 
pressure Hessian (Eq. (fTTj) ) in the i?*Q*-plane, as it was 
done in Ref. [||. Three main remarks can be made at 
this stage: (i) first, the pressure Hessian counteracts the 



development of the finite time singularity along the right 
tail of the Vieillefosse line implied by the RE term, (ii) 
probabilities are found to be very low in the dissipation 
production dominated region (i.e. R* > and above 
the Vieillefosse tail) meaning that pressure does not play 
there a significant role, and (iii) pressure redistributes 
the probabilities towards the enstrophy production dom- 
inated region (i.e. the flux is directed to the left, towards 
R* < 0). As far as the restricted Euler term is con- 
cerned, as is well known 0, Q , the deterministic Wre 
probability current pushes probabilities toward the right 
tail of the Vieillefosse line (data not shown). This re- 
sult helps create a picture of the time evolution of veloc- 
ity gradients along Lagrangian trajectories in stationary 
flows: The RE term "pushes" the probabilities towards 
the right tail towards and along the Vieillefosse line, while 
the pressure regularizes the implied finite time singularity 
and redistributes the probabilities towards the left part of 
the plane such that, in turn, the RE term can act again, 
etc. To that picture should be added the viscous diffusion 
effects, Wv pushes the probabilities toward vanishing R 
and Q not only along the Vieillefosse line but also ev- 
erywhere else, and stochastic forcing, such that the full 
probability current W (Eq. (jTOJ) ) is divergence free, in 
order to ensure stationary statistics (Eq. ©) 0. 

We display also in Fig. HJb) the amplitude of W p , us- 
ing logarithmic spacing of iso-probability lines. It can be 
seen that indeed no current is discernible in the right part 
of the plane. We remark that It would be interesting to 
quantify whether the effect of anisotropic pressure Hes- 
sian is orthogonal to the i?Q-plane in this region, when 
the RQ plane is extended into three-dimensions, as pro- 
posed and studied in Ref. [2(| ■ 



3. Incompressible Gaussian stochastic velocity field with 
K41 correlation structure 

Let us write a Gaussian homogeneous, isotropic and 
incompressible vectorial field u(x) [Il|,[22j], having a cor- 
relation structure consistent with K41 scalings. It reads, 
in d-dimensions, 



u e (x) 



m x - y) 



A dW(y) , (12) 



where dW(y) = (dWi(y), dW 2 (y), ...,dW d (y)) is a Gaus- 
sian vectorial white noise, and ipr, is a large-scale cut-off 
which involves the integral length scale L. The determin- 
istic kernel entering into Eq. (|12p is regularized over the 
small length scale e, namely |.| e = e * |.| (* stands for the 



convolution product), with a mollifier # € (x) 



(!) 



and J 6(x)dx =1. It is shown in Ref. [2l| that the ve- 
locity u e (x) has a well-defined limit when e — > 0, denoted 
by u(x), and such that (|u(x+£e)-u(x)|«) - C q (£/L)^ 3 
when i — > 0, with C q a constant independent on the vec- 
tor e. 

A Gaussian vectorial field, such as from Eq. (|T2"j) . is 
a poor representation of turbulence since it does not 
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FIG. 1: Probability current W P associated with the pressure Hessian -Eq. for the DNS velocity field and the three 

synthetic velocity fields -Eqs. (12)1. (TH) and JTSJ- in the 7T<2*-plane, where ft* = i?/ '(S^SV,) 3 / 2 and Q* = Q/(SijSi 3 ). The 
streamlines and vector plots of W P are shown in (a,c,e,g). The iso-probability contours of the magnitude of VV P are shown in 
(b,d,f,h). Contours are logarithmically spaced by factors of 10, starting at 1 for the contour closest to the origin. The thick 
lines represent the zero-discriminant (or Vieillefosse) line: ^-R 2 + Q 3 = 0. 



reproduce several important features such as a mean 
energy transfer towards small scales (i.e. the skew- 
ness phenomenon), the non-Gaussianity of velocity in- 
crements (i.e. the intcrmittency phenomenon) and the 
alignment of vorticity with the intermediate eigenvector 
of the strain rate tensor (JEHU- Nevertheless, it is useful 
to consider it in the analysis of the statistical quantities 
in which we are interested, such as the probability cur- 
rent W p associated with the pressure Hessian (Eq. (fTTIO . 
In particular, we look at which part can be attributed to 
Gaussian statistics and which part is really linked to tur- 
bulence. 

The Gaussian velocity field (Eq. (fT2)) ) is computed in a 
periodic box in d = 3 space dimensions, using N = 1024 3 
collocation points. The regularizing parameter e is cho- 
sen as e = 6dx, where the spatial resolution is dx = 1/N. 
For the mollificr 6 and the lar ge s cale cut-off ip, we take 
Gaussian functions. See Ref. 15| for further numerical 
details. The Pressure p is defined via the Poisson equa- 
tion Ap = — tr(A 2 ), where A is the (Gaussian) velocity 
gradient tensor. 

Figure [TJc) shows the vector plot and streamlines of 
the probability current W p obtained from the Gaussian 
velocity field -Eq. (fT2)K It can be seen that the pressure 
from the Gaussian field does not counteract directly the 
singularity along the direction of the right tail of the 
Vieillefosse line, as is the case in the DNS. Also, there is 
a significantly higher probability current in the right side 



region (i.e. R* > and above the Viellefosse tail) than in 
the DNS. The pressure obtained from the Gaussian field 
is only realistic in the enstrophy production dominated 
region (i.e. R* < 0), where the behavior shows indeed a 
trend to push the probability density towards this region. 
Also, the streamlines cross the Q = line meaning that 
the pressure field from the Gaussian velocity field does 
produce non-zero effect even when Q = 0. 

In Fig. [TJd), Wp-amplitude is shown. It can be seen 
that the amplitude iso- values are symmetric with respect 
to the R* = line. We arc thus led to the conclusion that 
a Gaussian velocity field, and its associated pressure field, 
do not make difference between dissipation production 
dominated regions (R* > 0) and enstrophy production 
dominated regions (R* < 0). 



Incompressible Multifractal stochastic velocity field with 
K062 statistics 



Based on the recent fluid deformation imposed by the 
Euler flow 0], and further heuristic introduction of the 
multifractal structure of turbulence as observed from ex- 
tensive empirical data (see e.g. [HI), Ref. 
the following 3D random vectorial field: 



15| proposed 



Ae Se(y) dW(y) , (13) 
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where S is a tensorial Gaussian log-correlated noise of 
the form 



-<T|<L 



(y - tr) ® [(y — cr) A rfW(tr)] 



,7/2 



[(y - cr) A dW(ff)] g> (y - <r) 



,7/2 



(14) 



with (g) denoting the tensorial product. The form of the 
symmetric matrix S e is inspired by the recent fluid defor- 
mation closure experienced by the fluid over short times 
[lH , and the exponent | has been selected such that the 
components of S are correlated logarithmically in space. 
A free parameter A enters this construction and governs 
the level of intermittency of the field. We will take in the 
sequel A 2 = 0.025 in order to be consistent with empirical 
findings 

Generation of the vectorial field u e (x) can be done 
accurately and efficiently in periodic boxes using up to 
1024 3 collocations points, in a similar way as done for 
the Gaussian velocity field (Eq. (fT2"j) ). The cost of the 
computation of the matrix exponential is the limiting nu- 
merical step. It is estimated at each point of space using 
a Pade approximant with scaling and squaring (see [15j 
for details). 

It has been shown numerically (Ref. [ll|) that the 
multifractal velocity field -Eq. (TT5)) - gives a realis- 
tic representation of instantaneous realizations of veloc- 
ity fields in fully developed turbulence in the inertial 
range, in regard to the following properties: (i) Longi- 
tudinal 8iu and transverse velocity increments are inter- 
mittent, A being the intermittency coefficient, (ii) the 
third-order moment ((5iu) 3 } is negative and proportional 
to the scale I. The fact that there is negative skew- 
ness S = {(Seu) 3 ) / {(Seu) 2 ) 3 / 2 means that u exhibits a 
non-vanishing mean energy transfer towards the small 
scales and (iii) vorticity gets preferentially aligned with 
the eigenvector of the strain-rate tensor corresponding to 
the intermediate eigenvalue. 

We show in Fig. QJe) the vector plot and streamlines of 
the probability current W p obtained from the multifrac- 
tal velocity field. It can be seen that, in a similar fashion 
to the Gaussian case, streamlines are roughly symmetric 
with respect to the R* = line. This is not consistent 
with DNS in the R* > region, but it is still realistic 
in the R* < region. The difference with the Gaussian 
field is the fact that now the joint density of R* and Q* 
is not symmetric with respect to the R* = line, show- 
ing thus the predominance of the regions for enstrophy- 
enstrophy production (upper-left quadrant, which in tur- 
bulent flows is correlated with vortex stretching) and 
dissipation-dissipation production (lower-right quadrant, 
connected with biaxial straining in turbulent flows). This 
is also the case for the probability current amplitude, as 
shown in Fig. [TJf). We can see therefore that both the 
Gaussian and Multifractal velocity fields do not repro- 
duce the void in probability in the R* > region as 
observed in DNS, but they do reproduce accurately the 



probability current evolutions in the R* < regions and 
the presence of probability flux at Q = 0. 



5. MTLM velocity field 

We consider a third case of a synthetic velocity field. 
The minimal turnover Lagrangian map (MTLM) veloc- 
ity field is obtained by distorting an initially random 
solenoidal vector field, uo(x), over a hierarchy of spa- 
tial scales {£ n — 2~ n L, n = 1, . . . , M}, where L is of the 
order of the integral scale, and the smallest scale, £m, 
is of the order of Kolmogorov scale. This generates the 
multiscale recursive sequence 



u„(x) = T[u„_i(x),4 



1, . 



(15) 



whose final step, Um(x), is the synthetic velocity field. 
Here, T[-] stands for the distortion operations applied. At 
each level n in the sequence the velocity is filtered at scale 
£ n and decomposed into low-pass and high-pass filtered 
parts: u< and u> respectively. The u < part is deformed 
by mapping the velocity vectors from their collocation 
points, x, to new positions that fluid particles moving at 
constant velocity in Lagrangian coordinates would reach: 
u<(X(t),t) = u<(x,0), withX(i) = x+£u<(x). The pa- 
rameter t is taken equal to the eddy-turnover time-scale 
corresponding to the spatial scale £ n , computed using 
standard Kolmogorov scaling. New velocity values at 
the collocation points are obtained by interpolation over 
nearby velocities that have come into a neighborhood of 
radius l n around x after the mapping. This deformed 
field u < is made solenoidal again by projection in Fourier 
space, and the amplitudes of its Fourier modes are scaled 
to conform to the target energy spectrum. Finally, u < 
is recombined with the u > part, which at this stage still 
remains as a Gaussian field. 

The next generation in the hierarchy will take the com- 
plete field u and will apply the same distortion opera- 
tions, now with the field decomposed at a smaller filtering 
scale. In this way, the effects are superposed and accu- 
mulated over a range of spatial scales. Further details 
and characteristics of these synthetic velocity fields can 
be found in Refs. [HI, EH- The present MTLM velocity 
field was generated in a periodic box, using 512 3 collo- 
cation points, with M = 6 generations in the hierarchy, 
and an energy spectrum corresponding to 1Z\ sw 250. 

The results for the probability current W p obtained 
from the MTLM velocity field (Eq. (|15p) are shown in 
Figs. [Hc-f). When compared with the DNS results (Figs. 
[IJa-b)), we can see close agreement of both magnitude 
and direction of W p . In particular, the MTLM velocity 
field reproduces, for pressure-related part of the prob- 
ability current, the void of probability in the R* > 
region. We can conclude that, of the three cases stud- 
ied, the MTLM velocity field gives the most realistic syn- 
thetic turbulence, as far as anisotropic pressure Hessian 
effects are concerned. But some small difference can be 
observed in the R* > region, close to the origin, where 
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FIG. 2: Conditional expectation (tr(P 2 )|Q*) of pressure Hes- 
sian norm on Q*\ DNS (solid), Gaussian (dot-dashed), Mul- 
tifractal (dashed) and MTLM (dotted) velocity fields. 



the MTLM fields seem to predict a circular motion that 
is not present in DNS. 



B. Mean pressure Hessian norm conditioned on Q 

As previously noticed in Ref. Q, current closures for 
the pressure Hessian 0, @, HH are proportional to the 
invariant Q. This is in particular the case for the closure 
provided in Ref. Q , namely 



— — Vtr(A 2 ) 



(16) 



where C Trj is the statistically stationary Cauchy-Green 
tensor at the Kolmogorov time scale T n 0, Q and 
tr(A 2 ) = — 2Q. Indeed, it is tempting to close the 
anisotropic part of P as a symmetric tensor proportional 
to Q, since the isotropic part is itself proportional to Q as 
seen on the Poisson equation Ap = 2Q. It would imply in 
particular that the probability current W p vanishes on 
the Q = line. This is not observed on DNS (see Figs. 
[TJa) and (b)). To quantify more precisely the behavior 
of P in the neighborhood of vanishing Q, we proposed in 
Ref. @ to estimate the average pressure Hessian (square) 
norm |P| 2 = tr(P 2 ) conditioned on the invariant Q. The 
corresponding conditional average (tr(P 2 )|<5) was shown 
for fields obtained by DNS as well as by applying the clo- 
sure (Eq. (HHP ). It was observed that for the DNS case, 
such conditional average at Q = does not vanish and 
furthermore it behaves nearly linearly with Q, whereas 
the closure Eq. (|16[) predicts a vanishing conditional av- 
erage for Q = and a quadratic behavior with Q. 



We display in Fig. [2]thc conditional average (tr(P 2 )|Q) 
as a function of Q for the four different velocity fields: 
the DNS, Gaussian (Eq. ([T2])). multifractal (Eq. (fT5|t ) 
and MTLM (Eq. (JTSJ) velocity fields. As previously 
observed, for the DNS case (solid line), the conditional 
average does not vanish at Q = and behaves linearly 
with Q in the neighborhood around Q = 0. For all the 
three remaining synthetic velocity fields, the conditional 
average does not vanish at Q — and hence they per- 
form better than the deterministic closure (Eq. [TB|) on 
this issue. For the Gaussian case however (dash-dotted 
line), the tails of the conditional average are not realis- 
tic, being far below the curves corresponding to the DNS 
case. Numerical simulations of the Gaussian fields at 
lower resolutions, i.e. TV = 256 3 or TV = 512 3 (data not 
shown), showed no difference with the TV = 1024 3 case. 
Interestingly, in this regard the multifractal field (dashed 
line) performed much better, exhibiting conditional aver- 
age tails very close to the DNS result. Some discrepancy 
is found for negative Qs where, the tail in the DNS case 
has a steeper slope. The MTLM velocity field (dotted 
line) also performs well against DNS data, although its 
tails are found to be quite symmetric, at odds with DNS. 

Overall, the behavior of the pressure Hessian obtained 
from the three synthetic velocity fields is reasonably sat- 
isfactory when compared against DNS data. The proba- 
bility current is well reproduced in the R < region and 
the conditional average shown in Fig. [5] does not vanish 
for vanishing Q. Only the MTLM velocity field can re- 
produce additionally the void in probability observed in 
DNS over the R* > region (explaining, or at least giv- 
ing an interpretation of the lack of action of the pressure 
Hessian in the dissipation production dominated region 
R* > remains, however, an open problem). Further- 
more, we have shown that synthetic velocity fields do 
predict the pressure Hessian square norm as being closer 
to linearly proportional to the invariant Q, rather than 
proportional to Q 2 as is the case in existing closures, in 
particular Eq. (|16|) . 

At this stage, one could reach the conclusion that some 
approximate surrogates of an actual turbulent field, even 
when obtained with the simplest Gaussian approxima- 
tion, contain a better prediction of the behavior of the 
pressure Hessian, in connection with its dependence on 
Q, than the deterministic closure given by Eq. p^|) . We 
will see in the following that the anisotropic part of the 
pressure Hessian can in fact be accurately closed by the 
local spatial variations of the invariant Q. 



III. LOCALITY OF THE PRESSURE HESSIAN 

We have seen in the first part of this work that a sim- 
ple Gaussian approximation, given by Eq. (|12p . or more 
sophisticated synthetic velocity fields, such as Eqs. (fT3|) 
and (|15|) . can reproduce the motion of the probability 
current W p associated with the pressure Hessian in the 
R* < region. Additionally, taking into account the 
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spatial distribution of the velocity field also leads to a 
non-vanishing conditional mean pressure Hessian norm 
for Q = 0. In this section we study to what degree spa- 
tial locality is important in determining these properties 
of the anisotropic part of the pressure Hessian. 

The exact expression @ for the pressure Hessian is 
very useful since it allows interpreting its isotropic part 
as being local, whereas the anisotropic part is governed 
by tr(A 2 ), or equivalently Q at different locations, i.e. it 
contains non-local contributions from the spatial varia- 
tions of Q. In this section, we will work with an equiv- 
alent form of Eq. ((2J that underlines the role played 
by the Hessian Q of the invariant Q. Indeed, taking 
two spatial derivatives of the Poisson equation that com- 
mute with the Laplacian, one obtains AP = 2Q, where 
Qij = a d r H . It is then easily seen that a similar rc- 

P and Q, 



J% 3 dxidx 

lation exists between the deviatoric parts of 



namely AP = 2Q , where the superscript d denotes the 
deviatoric part, i.e. for example P d = P — |tr(P)I, with 
I being the identity matrix. We finally reach a relation, 
equivalent to Eq. ([2]), between P d and Q d : 



P d (x) = - J- / ^—Q d (y)dy 
2?r J |x-y| 



(17) 



Relation (fl7|) is exact. In the following, we will truncate 
the integral present in Eq. (JTTJ) over a ball, centered at 
x, and of radius 77, namely 
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FIG. 3: In (a,b,c), we show the joint PDF of eigenvalues A p 
of the pressure Hessian and eigenvalues \ q of the proposed 
local expression (Eq. (|19|0 . Iso- probability lines correspond 
to e -4 , e -3 , e -2 , e _1 , 1. In (d) is shown the PDF of the cosine 
of the angle between the eigenvectors of P and those of the 
local closure (Eq. (|19[H : eigenvectors associated to the small- 
est (dashed), intermediate (solid) and biggest (dot-dashed) 
eigenvalues. 



- y \<v i x - y\ 



rQ d (y)dy 



(18) 



It is easily seen that from Eq. (|18p we recover Eq. p7|) 
by taking rj — > +00. We now make the strong assumption 
that T] is of order of the Kolmogorov length scale t\k ■ In 
this case, we can Taylor-expand the Hessian of Q at the 
position y around its value at the location x, take out 
Q d (x) from the integral, and perform the remaining in- 
tegration in spherical coordinates. We get an expression 
for the deviatoric part of the pressure Hessian 



d 2 p 

dxidx q 



( d 2 Q 

\dxidx-j 



(19) 



This now expresses the anisotropic part of the pressure 
Hessian in terms of local properties of the invariant Q, al- 
though the latter's spatial derivatives are needed. These 
derivatives are unknown in the closure and Lagrangian 
models of e.g. Q and so this does not constitute a prac- 
tical closure yet. To ensure that this expression yields 
the same norm as the true pressure Hessian, wc define 
the ball's radius according to 



(tr(P 2 )) 
(tr(Q 2 )) 



(20) 



The expression (fT!?|) and the choice of the length scale 
71 by Eq. (|2T)1) are consistent only if (i) eigenvalues and 



eigenvectors of P d and Q d are well correlated and (ii) r\ 
is indeed of the order of the Kolmogorov length scale tjk , 
since we assume that in the nei ghborhood of x, Q d (y) « 

Q d (x). 

Both Hessian tensors, for the pressure and for Q, are 
computed in Fourier space, for the periodic DNS flows. 
The expression (fH))) requires the computation of the sec- 
ond derivative of Q, which is already the square of a first 
spatial derivative. This is a reason for our use of the 
highly resolved DNS under study, in which dx/r/K ~ 1 
(see section |HX2]). Using r) K = {^/{e)) 1 ' 4 = 0.0123, we 
find that for the current DNS our choice of the length 
scale r] as Eq. (|20[) implies that rj = 1.83 r\K ■ Hence, 77 is 
of the order of r\K , as required. Let us stress that a well 
resolved DNS is necessary since some noise can be intro- 
duced by the computation of third order derivatives, lead- 
ing to a bad estimation of the parameter 77. More work is 
needed to clarify this point and to assess the dependence 
of 77 on resolution effects and Reynolds numbers. This 
could be done, for example, performing specifically de- 
signed DNS aimed at quantifying accurately high-order 
velocity derivatives, as was proposed in Rcfs. @, [H|. 

We display in Fig. [3] the joint probability densities of 
the three eigenvalues of the deviatoric part of the true 
pressure Hessian, (denoted by A p ) and the corresponding 
eigenvalues of the expression —r) 2 Q d (denoted by X q ). 
The smallest (Fig. 03a)) and largest eigenvalues (Fig. 
[3Jc)) are found to be well correlated with joint density 



8 




-1 1-1 o 1 




-1 1-2-1012 

R* -ti 2 Q> 2 

FIG. 4: (a) and (b): Probability current W P associated with 
pressure Hessian (Eq. for the DNS velocity field when 
using the local closure (|19l See caption of Fig. [TJ for further 
details, (c): Conditional expectation (tr(P 2 )|Q) of pressure 
Hessian norm on Q*\ DNS (solid) and closure (o), in a similar 
fashion than in Fig. [2] (d) Joint PDF of the components Pi 2 
and — rf AQ12 in a logarithmic scale. Isolines correspond to 
probability [10 -3 , 10" 2 , 10 _1 , 1, 10]. 



contours are elongated along the perfect correlation line 
(i.e. A p = \ q ), and correlation coefficients of p — 0.714 
and p = 0.703 respectively. 

As far as the intermediate eigenvalue is concerned, 
weaker correlation is found (p = 0.275), with isolines 
being close to circles. We can conclude that the smallest 
and largest eigenvalues are well correlated. In Fig. [3jd) 
we display the probability of the cosine of the angle be- 
tween the eigenvectors of the left and right terms of Eq. 
(flU)) . It is found that in all the cases, corresponding to 
the three different eigenvalues, the maximum of proba- 
bility is reached when the eigenvectors of P d and — rj 2 Q d 
are aligned. These results show then that P d and — rj 2 Q d 
are correlated both in amplitude and eigendirections. 

We have seen that, in tensorial structure, the true pres- 
sure Hessian P d and the local expression in terms of ve- 
locity gradients — rj 2 Q d are quite similar. We may now 
wonder if the local expression is able to reproduce the 
probability current associated to pressure effects as seen 
in Figs. QJa-b). To this purpose, we show in Figs. |4ja-b) 
the probability current W p obtained from DNS when the 
true pressure Hessian is replaced by the local expression 
(flU)) . We see in Fig. QTa) that indeed the local expres- 
sion reproduces the counteractive action of the pressure 
along the right tail of the Vieillefosse line. Furthermore, 
a void in probability is found in dissipation production 
dominated regions (i.e. R* > 0), as can be clearly seen in 



Fig. Qlb). Finally, the local expression also reproduces 
some of the probability redistribution in the enstrophy 
production dominated region (R* < 0), as observed in 
DNS (Fig. nja-b)), although the direction of the prob- 
ability flux streamlines for R* < is seen to be more 
vertical than the left-wards directions seen in the DNS. 
Also, the streamlines are found more curved for the local 
expression than for P and the probability current am- 
plitude |Wp| is found to decrease faster at high values 
of R and Q than in the DNS case. This could be due 
to limitations of the localized expression in reproducing 
very high turbulent fluctuations 

The overall behavior of the local expression in the RQ- 
plane is on the whole quite satisfactory when compared 
against DNS. Some differences appear: (i) the stream- 
lines of the probability current are found more curved 
for the local closure than for P, very much in opposition 
to the streamlines imposed by the RE approximation (see 
Ref. [1]), and (ii) the probability current amplitude |W P | 
is found to decrease faster at high values of R and Q than 
in the DNS case. This could be due to limitations of the 
localized expression in reproducing very high turbulent 
fluctuations. On the whole, however, the trends pro- 
vided by the local expression Eq. (fT5)) agree quite well 
compared to the DNS results in terms of the prbability 
fluxes in the RQ plane. 

We also show in Fig. |4jc) the conditional expectation 
of the pressure Hessian square norm, based on Q. For 
the sake of clarity, we show again the conditional aver- 
age as obtained in DNS (solid line) (it was already shown 
in Fig. [2|). We use open symbols (o) to show the con- 
ditional expectation obtained from DNS using the local 
approximation Q19p . It can been seen that the conditional 
expectation obtained from the local approximation is in 
very good agreement with the DNS, with some discrepan- 
cies appearing for negative Q* . Interestingly, we see that 
the conditional expectation does not vanish for Q = 0. 
This is consistent with former remarks made about the 
requirement that a realistic model of P cannot be simply 
proportional to Q. 

Finally, to quantify the agreement between individual 
tensor elements, we show in Fig. 0|d) the joint PDF of 
the component P±2 and the component — ?7 2 Qi2- This 
plot demonstrates the good level of correlation between 
these two tensors since the joint PDF is clearly elongated 
along the perfect correlation line (i.e. P12 = —T] 2 Qi2)- 
The corresponding correlation coefficient is p = 0.55. 
Thus, also this statistical test confirms the good agree- 
ment between the local approximation and the true pres- 
sure Hessian. 



IV. CONCLUSIONS AND PERSPECTIVES 

This article focuses on the statistical nature of the pres- 
sure Hessian that governs much of the Lagrangian dy- 
namics of the velocity gradient tensor in turbulence. In 
a first part, we have seen that synthetic velocity fields re- 
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produce many properties of the pressure Hessian as they 
are seen in DNS flows, such as the non-trivial behavior of 
the probability current, and the conditional expectation 
of the pressure Hessian norm on the invariant Q. Even 
the simplest Gaussian approximation for the velocity field 
(Eq. (fT2j) ) can represent some non-trivial behaviors of the 
pressure that could not be predicted by the closures in 
terms of A 0,0, EH- Based on this observation and on an 
exact field description of the pressure Hessian by means 
of nonlocal integrals (Eqs. © and (fTTf b we formulate 
the hypothesis that considering only the integration over 
a ball of radius t]k, and neglecting other contributions, 
the deviatoric part of P could be expressed in terms of 
local properties of the velocity gradient tensor, but in 
terms of higher-order derivatives. Specifically, the spa- 
tially local approximation is not expressed in terms of A 
but in terms of second-order derivatives of Q. This ap- 
proximation was found to be highly correlated with the 
true pressure Hessian P when compared in DNS compu- 
tations. These findings show that the main contribution 
to P(x) is contained in the local neighborhood around 
position x, in a ball centered at x and of radius of the 
order of tjk ■ This raises the hope that local closures in- 
volving a finite set of ordinary differential equations may 
still be possible, for studying the Lagrangian dynamics of 
the velocity gradient tensor. To that end, it is still nec- 
essary to express the Hessian of Q in terms of the local 
values of A. Only then would we have a full closure. 

Let us finally remark that if a tractable transport equa- 
tion for the pressure Hessian is difficult to get, the La- 
grangian derivatives of p and P can be related. We get 
the following transport equation for the pressure Hessian 



d°P 

dt 



.^vv. + vvf. ( 21 , 



tive or Oldroyd derivative, that relates the rate of change 
written in the coordinate system rotating and stretching 
with the fluid (see for example Ref. [H]), be. 



d°P dP Tt , 
—— = — + A T P 

dt dt 



PA 



(22) 



Then, if we neglect the right-hand side of Eq. (|2"Tj) in 
the rate of change of the pressure Hessian, i.e. if we 
assume that the Oldroyd derivative (f2"2"j) vanishes, and we 
apply the recent fluid deformation (RFD) approximation 
(assuming A independent on time), we get that for an 
early time 



P(r) 



-rA 



P(0)e 



-tA 



(23) 



where d /dt stands for the upper convected time deriva- 



This remark justifies the use of matrix exponentials as 
closures of the pressure Hessian Q, as it was also noted 
for the subgrid-scale stress tensor [27j ■ If in addition 
we start from an initial isotropic pressure Hessian, and 
include the Poisson equation, we recover exactly the clo- 
sure ([16]) proposed in Ref. This represents useful 
perspectives for future investigations on the Lagrangian 
dynamics of the velocity gradient tensor. 
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